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ABSTRACT 

A scale-invariant model of statistical mechanics is applied to present invariant forms of 
mass, energy, linear, and angular momentum conservation equations in reactive fields. The 
resulting conservation equations at molecular-dynamic scale are solved by the method of large 
activation energy asymptotics to describe the hydro-thermo-diffusive structure of laminar 
premixed flames. The predicted temperature and velocity profiles are in agreement with the 
observations. Also, with realistic physico-chemical properties and chemical-kinetic parameters 
for a single-step overall combustion of stoichiometric methane-air premixed flame, the laminar 
flame propagation velocity of 42. 1 cm/s is calculated in agreement with the experimental value. 

1. INTRODUCTION 

The universality of turbulent phenomena from stochastic quantum fields to classical 
hydrodynamic fields resulted in recent introduction of a scale-invariant model of 
statistical mechanics and its application to the field of thermodynamics [4]. The 
implications of the model to the study of transport phenomena and invariant forms of 
conservation equations have also been addressed [5]. In the present study, the invariant 
forms of the conservation equations are described and the results are employed to 
introduce a modified hydro-thermo-diffusive theory of laminar premixed flames. 

2. A SCALE-INVARIANT MODEL OF STATISTICAL MECHANICS 

Following the classical methods [1-3], the invariant definitions of the density pp, and 
the velocity of atom Up, element Vp, and system w^ at the scale {3 are given as [4] 


pp=npmp=mp| fpdup , > v p=Pp m p| upfpdup 


w rVi 


( 2 . 1 ) 


The scale-invariant model of statistical mechanics for equilibrium fields of . . . eddy-, 
cluster-, molecular-, atomic-dynamics ... at the scale (3 = e, c, m, a, and the 
corresponding non-equilibrium laminar flow fields are schematically shown in Fig.l. 
Each statistical field, described by a distribution function fp( Up) = fp(rp, Up, tp) drpdup. 
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defines a "system" that is composed of an ensemble of "elements", each element is 
composed of an ensemble of small particles viewed as point-mass "atoms". The 
element (system) of the smaller scale (P) becomes the atom (element) of the larger 
scale (P+1). The three length scales associated with the atom, the element, and the 
system at any scale (3 are (lp = 2^, A.p, Lp = Xp +1 ) where A.p = <l 2 p> 1/2 is mean-free - 
path of the atoms [5]. 
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Fig.l Hierarchy of statistical fields for equilibrium eddy-, cluster-, and 
molecular-dynamic scales and the associated laminar flow fields. 

The invariant definitions of the peculiar and the diffusion velocities have been 
introduced as [4] 

v 'e = u H- v fi . v p = v p- W P ( 2 . 2 ) 

such that 

Vp = V'p +1 • (2.3) 

The above definitions are applied to introduce the invariant definitions of equilibrium 
and non-equilibrium thermodynamic translational temperature and pressure as [4] 
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3kTp = mp<up 2 > 

, Pp = pp<up 2 >/3 

(2.4) 

and 



3kTp= mp<V' 2 p> 

, Fp = npmp<V' 2 p>/3 

(2-5) 

leading to the corresponding 

invariant ideal "gas " laws [4] 


PpV = NpkTp and 

F pV — NpkT p 

(2.6) 

3. SCALE-INVARIANT FORM OF THE CONSERVATION EQUATIONS 

FOR 


CHEMICALLY REACTIVE FIELDS 

Following the classical methods [1-3], the scale-invariant forms of mass, thermal 
energy, linear and angular momentum conservation equations [5] at scale (3 are given as 


dpp 

dt 

+ v. (ppVp) = Qp 

(3.1) 

d£p 

dt 

+ V. (epVp) = 0 

(3-2) 

dt 

+ V. (ppVp) = 0 

(3.3) 

3itp 

at 

+ V. (jtpVp) = 0 

(3.4) 


where = pphp, pp = PpVp, and ftp = PpCOp are the volumetric density of thermal energy , 
linear and angular momentum of the field, respectively and to p = Vxvp i s the vorticity. 
Also, £2p is the chemical reaction rate, and hp is the absolute enthalpy [5]. 

The local velocity Vp in (3.1)-(3.4) is expressed as the sum of convective Wp = <Vp> 
and diffusive velocities [5] 


V e = w p + v pg 

5 

^pg - - Dp v i n (pp) 

(3.5a) 

V p=Wp+ Vptg 

? 

V ptg =-apVln(sp) 

(3.5b) 

V p=Wp+Vpb g 


v phg ~ ~ v p V ln(Pp) 

(3.5c) 

Vp=Wp+ Vprhg 

y 

Vprhg = ~ v p ^ ln(3Tp) 

(3.5d) 


where (Vpg, Vg tg , Vph g , Vprhg) are respectively the diffusive, the thermo-diffusive, the 
linear hydro-diffusive, and the angular hydro-diffusive velocities. For unity Schmidt and 
Prandtl numbers Scp= Prp = Vp/Dp = Vp/otp = 1, one may express 



* prh g ~ Vpg + Vprh 

that involve the thermal V Bt) the linear (translational) hydrodynamic and the angular 


(rotational) hydrodynamic V^-h diffusion velocities defined as [5] 

Vpt =-a p Vln(h p ) (3.6a) 

Vpi = - Vp V In(v p ) (3.6b) 

V j)th = -VpVln(cOp) (3.6c) 

Since for an ideal gas hp = CppTp, when Cpp is constant and T = Tp, Eq.(3.6a) reduces to 
the Fourier law of heat conduction 

q=phpVp t = - KpVT (3.7) 

where Kp and cu = Kp((pp c pp) are the thermal conductivity and diffusivity. Similarly, 
(3.6b) may be iaentified as the shear stress associated with diffusional flux of linear 
momentum and expressed by the generalized Newton law of viscosity [5] 

p = PpVjpVijph = - Up Vv p = - pp avjp/axi (3.8) 

Finally, (3.6c) may be identified as the shear stress induced by diffusional flux of angular 
momentum (torsional stress) and expressed as 

Fjr p = Pp^jpVijprh = — Pp = ~ P'P ^jp/dXi (3.9) 

Substitutions from (3.5a)-(3.5d) into (3.1)-(3.4), neglecting cross-diffusion terms 
and assuming constant transport coefficients with Scp = Prp= 1, result in 

+ wp . Vp p - DpV 2 pp - Qp (3.10) 


hp 


^ + wp . Vpp - DpV 2 pp 
0t 


+ PP 


dh| 

It 


£■ + wp . Vhp - apV 2 hp 


= 0 


(3.11) 


Vp 


% + wp . Vpp - DpV 2 pp 
at 


cop + wp . Vpp - DpV 2 pp 
8t 


+ pp 


+ Pp 


dv 


at 


a© 


+ wp . Vvp - vpV 2 vp 


0 


[ at 


£• + wp . Vcop-t- top . Vwp- vpV 2 cop 


(3.12) 


= 0(3.13) 


The above forms of the conservation equations perhaps help to better reveal the 
coupling between the gravitational versus the inertial contributions to total energy and 
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momentum densities of the field. Except for possible externally imposed sources, Ep, pp 
and Kp have no internal sources as reflected in (3.2)-(3.4). However, in the presence of 
chemical reactions the loss of gravitational mass could result in the production of inertial 
thermal energy, or linear and angular momenta. For example, the first and the second 
parts of (3.11) respectively correspond to the gravitational and the thermal 
contributions to the total energy density of the field. For instance, the loss of 
gravitational mass induced by chemical reaction in the body of a person results in the 
generation of thermal energy (heat) in this persons body. Similarly, the first and the 
second parts of (3.12) respectively correspond to the gravitational and the inertial 
contributions to the total linear momentum density of the field. Now, one considers a 
stationary person with no initial linear momentum that suddenly starts to run, thus 
producing substantial linear momentum without the action of any external forces. In 
this case, there is no violation of the conservation of momentum, but rather because of 
chemical reactions in the body of such a person, the first part of (3.12) changes thus 
leading to a compensating change in the second part. Finally, the first and the second 
parts of (3.13) respectively correspond to the gravitational and the inertial 
contributions to the total angular momentum density of the field. For example, (3.13) 
may be used to describe the change of angular velocity of a ballet dancer. Here, the 
loss of mass by chemical reactions in the body of a spinning dancer that pulls the arms 
inwards, thus doing work against centrifugal forces, leads to an increase in the dancer's 
angular momentum. Because of the large value of the velocity of light c in E = me 2 , the 
actual loss of gravitational mass in the above examples will be exceedingly small. 

Substitutions from (3.10) into (3.11)-(3.13) result in scale-invariant forms of 


conservation equations [5] 


^ + wp ■ Vpp - DpV 2 pp = Qp 

at 

(3.14) 

^12. + wp . VTp - apV^p = - hpQp/(ppc p p) 
at 

(3.15) 

— + wp . Vvp - vpV 2 vp = - vpQp/pp 
dt 

(3.16) 

<9(0 g 2 

+ Wp . Vcop - vpV a) p = - (Op. Vwp - oopQp/pp 

(3.17) 


Equation (3.17) is the modified form of the Helmholtz vorticity equation for chemically 
reactive flow fields. The last two terms of (3.17) respectively correspond to vorticity 
generation by vortex-stretching and chemical reactions. Also, equation (3.16) is the 
scale-invariant equation of motion in reactive fields [5] that includes the reaction term 
(-vpQp/pp) representing generation Qp < 0 (annihilation Qp > 0) of linear momentum 
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accompanied by release (absorption) of thermal energy associated with exothermic 
(endothermic) chemical reactions. It is known that as flames propagate, they convert 
stationary reactants to moving combustion products because of thermal expansion. 
Another important feature of the modified equation of motion (3.16) is that it involves a 
convective velocity wp that is different from the local fluid velocity vp. Consequently, 
when the convective velocity vanishes wp = 0, equation (3.16) reduces to the diffusion 
equation similar to mass and heat conservation equations (3.14)-(3.15). Because the 
convective velocity wp is not locally-defined it cannot occur in differential form within 
the conservation equations [5]. This is because one cannot differentiate a function that 
is not locally, i.e. differentially, defined. To determine wp, one needs to go to the next 
higher scale 0+1) where Wp= vp+i becomes a local velocity. However, at this new 
scale one encounters yet another convective velocity wp + i which is not known, 
requiring consideration of the higher scale 0+2). This unending chain constitutes the 
closure problem of the statistical theory of turbulence discussed earlier [5]. 

By summation of (3.1)-(3.4) over O) one can arrive at the conservation equations 
at the next higher scale of 0+1). By such procedure, one can move from molecular- 
dynamic to cluster-dynamic scale or from cluster-dynamic to eddy-dynamic scale 
within the cascade of embedded statistical fields (Fig.l). The summation of Eq.(3.1) is 
simple since 

2ppp=pp+i (3.18) 

and 

2p ppvp = 2p [pp+iYpVp] = pp+i 2p [Ypvp] = pp+ivp+i (3.19) 

For Eq.(3.3), the summation of the first term is identical to that shown in (3.19) 

2ppp = 2pppvp= 2p [pp +1 YpVp] = pp + i 2p [Ypvp] =pp+iVp+i= pp+i (3.20) 

To treat the summation of the second term of (3.3), one starts with the relation based on 

(2.1)-(2.2) 

vp = wp + Vp = vp + i + V'p+i (3.21) 

Multiplying (3.21) by (Y p + i ppvp) and summing over (P) and (P+1) leads to 

2p + i 2p [Yp+ippvpvp] =2p+i2p [Yp+ippvpvp + i] + Ip+i 2p [Yp +1 ppVpV'p + i] 


or 


2p [ppvpvp] = 2p + i2p [ppVpYp + iVp +1 ] + 2p + i 2p [ppVpYp + iV'p + i] = 

= [ppvpwp + i ] + 2p [ppVpVp + i] = 

= Pp + ivp + iwp + i + pp+iVp +1 Vp + i = 

= Pp+i v p+i( w p+i + Vp + i) = pp+iVp + iVp +1 = pp+ivp +1 (3.22) 



7 


where Yp is mass fraction and use was made of the relation vp + i = wp+i + Vp+i from 
(2.2) in the last step. The summation of the energy (3.2) and vorticity (3.4) equations 
follow procedures similar to those used above in (3.21)-(3.22). 

4. CONNECTION BETWEEN THE MODIFIED FORM OF EQUATION OF 
MOTION AND THE NAVIER-STOKES EQUATION 

The original form of the Navier-Stokes equation with constant coefficients is given 
as [1, 2] 


p — + pv . Vv = - VP + pV 2 v + i pV(V.v) 
dt 3 


(4.1) 


when the coefficient of bulk viscosity is considered to be negligible [6]. Since 
thermodynamic pressure P t is an isotropic scalar, P in (4.1) is not P t . Rather, the 
pressure P is generally identified as the mechanical pressure that is defined in terms of 
the total stress tensor Tij = - P t Sjj + x,j as [6] 


P m =-(l/3)Tii = P t - (1/3) xji (4.2) 

The normal viscous stress is given by (3.8) as (1/3) Xii = (1/3) pviVii = - (1/3) pV.v and 
since VP t ~ 0 because of isotropic nature of P t , the gradient of (4.2) becomes 


VP = VP m = V i p(V.v ) = i pV ( V.v ) (4.3) 

Substituting from (4.3) in (4.1), the Navier-Stokes equation assumes the form 

— + v . Vv - vV 2 v = 0 (4.4) 


that is almost identical to the modified equation of motion (3.16) with Qp = 0 except that 
in the latter the convective velocity Wp is different from the local velocity Vp. However, 
because (4.4) includes a diffusion term and the velocities Wp and Vp are related by Vp = 

Wp + Vp, it is clear that (4.4) should in fact be written as (3.16). 

An example of exact solution of the modified equation of motion (3.16) was recently 
introduced [7] for the classical Blasius problem [2] of laminar flow over a flat plate. For 
this steady problem Eq.(3.16) in the boundary layer, with w'y = 0 and Q = 0, reduces to 


w 


dv'x 

dx' 


dy' 2 

v 'x = 0 

v'x = w'o 


(4.5) 


y' = 0 

y' = oo 


= v 


(4.5a) 

(4.5b) 
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where w' Q is the constant free-stream velocity outside of the boundary layer and (x', y’) 
are the coordinates along and normal to the wall, respectively. The local velocity v' x 
varies from v' x = 0 at the wall y' = 0 to v' x = w' G at the edge of the boundary layer at all 
axial positions. Therefore, the convective velocity w' x = <v' x >, i.e. the mean value of the 
local velocity v' x inside the boundary layer, will have the constant value of w' x = w'q/2. 
Substituting for this convective velocity in (4.5) one obtains [7] 


211 % 

o 

II 

* (S 
+ 

(4.6) 

T] = 0 

v x = 0 

(4.7a) 

r\-co 

v x = 1 

(4.7b) 


in terms of the similarity variable r\ = y7(2x' 05 ), where v x = v' x /w’ 0 . The solution of 
(4.6)-(4.7) is the predicted velocity profile v x = erf q that is in excellent agreement with 
the experimental observations of Nikuradze [2], 

The solution of (3.16) for the classical problem of Hagen- Poiseuille flow [2] in 
circular tubes has also been investigated [7]. It was found that the geometry of the 
predicted velocity profile involving Bessel function was quite similar to the classical 
parabolic profile and hence in agreement with the experimental observations. Finally, 
the exact solution of the modified Helmholtz vorticity equation (3.17) was recently 
reported [8] for the steady problem of non-reactive flow within a stationary liquid droplet 
that is located at the stagnation-point between two axisymmetric counter-flowing 
gaseous streams. It was found that such a spherical flow could be expressed by the 
stream function [8] 


*F = - r 2 z (1 - r 2 - z 2 ) 


(4.8) 


representing two rings vortices, rather than a single spherical Hill vortex for the 
classical problem of a droplet in a uniform stream [6], that are located above and below 
the stagnation plane. Therefore, the preliminary investigations discussed above show 
that the modified equation of motion (3.16) does indeed lead to realistic solutions in 
agreement with experimental observations for these classical problems for which exact 
solutions of the Navier-Stokes equation are available. 

5. MODIFIED HYDRO-THERMO-DIFFUSIVE THEORY OF LAMINAR 
FLAMES 

Theory of laminar flames is the most fundamental problem of combustion science 
and subject of many classical [9-18] as well as more recent [3, 19-24] studies. For one- 
dimensional propagation of a planar flame one introduces the dimensionless parameters 

A * [v F W F Ba/(pv o)] e- fa (5.1) 


6 = (T - T u )/(T b - T u ) , 


y = Y F /Y Fu 



The adiabatic flame temperature Tb, the Zeldovich number p, and the coefficient of 
thermal expansion x are 

T b = T u + QY Fu /(v f Wfc p ) , P = E(T b - T u )/RT b 2 , X = Ob ~ T u )/T b (5.2) 

and one assumes that p »1. Also, Q, E, and R are heat release per mole of fuel, 
activation energy, and universal gas constant. Prandtl, Schmidt, and Lewis numbers 
are assumed to be unity v = a = D, such that outside of reaction zone where A = 0, the 
0, y, and v fields will be similar under identical boundary conditions. Equations (3.14)- 
(3.16) for molecular-dynamic scale with the dimensionless coordinate, time, and 


velocity 

defined as x = 

x7(a/v' 0 ), t = 

t'/(a/v'2 0 ), and v = v7v' 0 become 



ay 

w — = 

d 2 y 

-Ay eP( 0 ~ 

A 6(x) 

(5.3) 

at 

dx 

dx 2 



d A + 

ae 

w = 

a 2 e 

+ A y eP (e ~ 

V 6(x) 

(5.4) 

at 

dx 

dx 2 



+ 

dv 

w — = 

a 2 v 

+ v A y eP< 6 

- A 5(x) 

(5.5) 

at 

dx 

dx 2 




where v' 0 is the flame propagation speed. 

5.1 Far-Field Convective Coordinate 

For laminar flames propagating in quiescent reactive fields, there is no forced 
convection w = 0, and (5.3)-(5.5) reduce to non-homogeneous diffusion equations with 
nonlinear sources. Because of thermal expansions in the flame, the stationary cold 
reactants are converted to moving hot combustion products resulting in a velocity jump 
across the flame sheet as schematically shown in Fig.2a. When viewed from the 
perspective of the physical or far-field coordinate x\ the flame appears as a 
mathematical surface of discontinuity without any spatio-temporal structures. 


X '=-00 X 1 »0 x' = 00 z' = - CO t = 0 z'=°° 
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Fig.2a A propagating laminar flame Fig.2b A stationary laminar flame 

viewed from far-field coordinate x\ viewed from far-field coordinate z'. 
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Outside of the thin flame zone chemical reactions will be frozen A = 0, and the governing 
equations (5.3)-(5.5) become 


df 

dt' 


= v 


d 2 f 


m 


dx 


1 2 


f = y, 0, v 


(5.6) 


To an observer in the coordinate x' only the mean velocity w'f = - v' 0 - v't/2 will be 
detectable and this velocity will be without any spatial structure. If one introduces a 
coordinate system that moves with the flame z’ = x' - w'f t', one obtains from (5.6) 


Wf 


= . v dH 


dz' 


m 


dz 


i2 


f = y, 0, v 


(5.7) 


and the flame becomes stationary with velocity jump (~v' 0 , -v' 0 -v'b) across the flame 
front (z'= 0+, z'= 0 -) as shown in Fig.2b . However, the hydro-thermo-diffusive structure 
of the flame will be hidden from this observer in the far-field coordinate z'. 

5.2 Outer Convective-Diffusive Coordinate 

In terms of the convective-diffusive stretched coordinate x = x’/It, where It is the 
flame thermal thickness \j= a/v' 0) the hydro-thermo-diffusive structure of the flame 
becomes visible as shown in Fig.3a. However, the convective-diffusive coordinate is 
being applied to the colder regions outside of the thin reaction zone such that chemical 
reactions remain frozen A = 0 because (3 »1. With linear approximation, the values of 
velocities (v, w, V) within the flame structure at the locations x = (- 1/2, 0, 1/2) will be 
v = ( - v b , - Vb/2, 0), w = (- 3vb/4, - Vb/2, - v\J4\ and V = (- Vb/4, 0, Vb/4) as shown in 
Fig.3a, when Vb= v'b/v' 0 . Hence, the flame hydrodynamic structure involves a variable 
convective velocity w, and the relation v = w + V is satisfied, while the sign of V as well 
as that of curvature of the velocity profile change across x = 0. 



Fig.3a Hydrodynamic structure of a 
propagating laminar flame viewed 
from thermo-diffusive coordinate x. 



Fig.3b Hydrodynamic structure of a 
stationary laminar flame viewed from 
thermo-diffusive coordinate z. 
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The velocity of the burned gas v'b will be related to the laminar flame propagation 
velocity v' Q by the mass conservation across the flame p 0 v' 0 = Pb(v'o+ v'b) that is valid to 
a very high degree of accuracy. It is important to emphasize however that mass is not 
strictly conserved across flame fronts and in fact an exceedingly small amount of 
gravitational mass (rest energy) of the reactants will be converted to other forms of 
energy (inertial energy) such as thermal energy of products of combustion and photons 
as discussed in sec. 3. 

For the steady problem, the convective-diffusive coordinate that moves with the 
flame is z = x - w t where w = wYv' 0 , z = z71 T , and t = t7(l T /v'o). One notes however 
that in terms of the stretched coordinate z, the velocity w = - 1 + Vb(z - l)/2 is no 
longer constant but has spatial structure as shown in Figs.3a-3b. One further notes that 
the average value of w at the flame center z = 0 is Wf = - 1- vt/2 that is exactly the 
mean flame speed as seen by the far-field coordinate z' discussed in sec.5.1 above. 
Furthermore, for an observer that moves with the flame front, this velocity profile 
w = <v> = - 1 + Vb(z - l)/2 remains invariant such that Eq.(5.7) becomes 


[1 - v b (z - l)/2] * = f = y, 9, V 

dz dz 2 

Z -* 00 y— 1 = 9 = V + 1 = 0 

z - <x> y=9-l=v + l + Vb = 0 

Introducing the new variable fvb" t, = 1 - Vb(z - l)/2 into (5.8)-(5.10) gives 

2 S# + ^ = 0 f= y,6,v 

dt 2 

£ — » - oo y-l=9=V+l=0 

^ oo y = 0- l= v+ l + Vb = O 

with the solutions 


(5-8) 

(5.9) 

(5.10) 


(5.11) 

(5.12) 

(5.13) 


0 = 1 - y = - (v +l)/v b = (1 + erf £)/2 (5.14) 

According to (5.14), the steady flame structure is given by a traveling error function 
schematically shown in Fig.4b. This result, while being in harmony, is fundamentally 
different from that discussed in an earlier study [25], and perhaps provides a clearer 
description of the steady flame structure. The predicted geometry of the temperature 
profile involving error function (5.14) is in close agreement with the experimental 
observations [26-28]. In addition, the predicted profiles are consistent with the 
measured temperature profiles of counterflow premixed flames [29] in the limit of 
vanishing rates of stretch. On the other hand, the temperature profile within the flame 
structure according to the classical theories of laminar flame is given by an exponential 
function as schematically shown in Fig.4b. More experimental measurements of 
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temperature profiles in unsteady propagating laminar flames are needed in order to 
further test the validity of the modified versus the classical theory. The slope of the 
temperature profile of the outer convective-diffusive zone at the position of the reaction 
zone £i (Fig. 4a) is obtained from (5.14) as 


|de] _ siSt (5.15) 

that will be matched with the solution within the reaction zone to be described next. 



Fig.4a Flame structure according to 
the modified theory of laminar flame. 



Fig.4b Flame structure according to 
the classical theories of laminar flame. 


5.3 Inner Reactive-Diffusive Coordinate 

The analysis of the thin reaction zone follows the classical methods [3, 19-24] and 
involves the stretched coordinate 

S = 0(£i-o (5.16) 

along with the temperature and concentration expansions 

8=l + 0i/p + ... and y = 0 + Yi/p + ... (5.17) 

that are introduced into Eqs.(5.3)-(5.4) to obtain, to the first order in e = l/j3 « 1 

fXi = -4/L. Y, e®' (5.18) 

d§ 2 v b p 2 

d 2 ®i — _ 4 A y 1 e 01 
v b p 2 


(5.19) 
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From the coupling of Eqs.(5.18)-(5.19) and the boundary conditions at 
obtains ©1 + Yi = 0 such that Eq.(5.19) becomes 


d20 l = 4 A 0 e 0. (5.20) 

d| 2 v b p 2 

that may also be expressed as 

d /d0 i\ 2 _ 8 A 0 e e, (5.21) 

d0i v d^ / v b p 2 

The first integral of the above equation, and matching of the slopes of the temperature 
profiles on either side of the reaction zone with the outer solutions in Eq.(5.15) results in 


exp (- 2 If) = [8jt v F W F Ba/(pv b v o P 2 )] e" p/x (5.22) 

that relates the reaction zone position ^ and hence the ignition temperature ©i, to the 
flame propagation velocity v' Q . The parameter B in (5.1) is related to the actual 
preexponential factor B' in the law of mass action [3] under Arrhenius kinetics by 

B = B’ppo/(W F Wo) (5.23) 

Also, the mass balance across the flame front p u v' 0 = p b (v' 0 + v' b ) leads to 

V' b = V' 0 (pu - p b )/pb ( 5 - 24 ) 

By substitution from (5.23)-(5.24) into (5.22), one obtains the analytic expression 


v-2 _ v F p b po a B' g - p/ x ^ (5 25) 

(pu-Pb)WoP 2 

for calculation of laminar flame propagation velocity. 

For single-step overall combustion of stoichiometric premixed methane-air flame at 
the flame temperature of 2100 K the relevant physico-chemical properties are v F = 1, 
po = 1 .38 kg/m 3 , Wo = 32, a ~ 6x 10- 4 m 2 /s (thermal diffusivity of air at the flame 
temperature of 2100 K), E ~ 46 kcal/mole , B' ~ 4.33xl0 7 m 3 /kmol-s [30], x, ~ 0.86, 

P ~ 10, and the ignition temperature of 0i = 0.99 that by (5.13) gives the reaction zone 
position = 1.75. Also, from the temperatures of reactants 300 K and combustion 
products 2100 K and the ideal gas law under constant pressure, one obtains the density 
ratio p u /p b = 7. With these realistic values of the physico-chemical properties, the value 
of flame propagation velocity calculated from (5.25) is about v' Q = 42.1 cm/s in close 
agreement with the experimentally observed value [27, 28, 31-32]. Although this level 
of agreement between the theory and experiments is considered to be encouraging, it 
should be viewed with caution because of the well-known uncertainties in the overall 
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chemical-kinetic parameters (E, B’)- The value of about v' Q = 42 cm/s has also been 
obtained in a number of numerical investigations using complex multi-step kinetic 
models [33-36]. 

6. CONCLUDING REMARKS 

A scale invariant model of statistical mechanics was applied to present invariant 
forms of mass, energy, linear, and angular momentum conservation equations in 
chemically-reactive flow fields. The summation procedures for relating adjacent 
families within the hierarchy of statistical fields was described. Also, the coupling 
between the gravitational versus the inertial contributions to the energy-momentum 
density of the field was discussed. The connection between the modified equation of 
motion and the classical Navier-Stokes equation was established. The exact solution of 
the modified equation of motion for the classical problem of Blasius for laminar flow over 
a flat plate was presented. Also, a modified form of the Helmholtz vorticity equation 
was presented with a source of vorticity due to chemical reactions. The conservation 
equations at the molecular-dynamic scale were then applied to present a modified 
hydro-thermo-diffusive theory of laminar flames. The predicted flame structure was 
found to be in agreement with experimental observations as well as numerical 
calculations. With realistic physico-chemical properties for one-step combustion of 
stoichiometric methane-air premixed flames, the flame propagation velocity of 42.1 cm/s 
was calculated in accordance with experimental observations. 
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